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Abstract 

We study the stability of uniformly moving membrane-like objects in seven 
dimensional Anti-de Sitter space. This is approached by a linear perturbation 
analysis and a search for growing modes. We examine both analytic and 
numerical configurations previously found in 
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1 Introduction 

The paper [T] recently explored the properties of defects moving through 
Anti-de Sitter space. This was a generalization of previous works, [2 [3], 
which examined the important case of a moving string in AdSa . The focus of 
[1] was on the case of a membrane in AdSy for reasons such as ease in analytics 
and relevance from a quantum gravity (M-theory) point of view. The major 
use of such explorations is in the context of the AdS/CFT correspondence [1]. 
This duality allows properties of the defects, as calculated in the low energy 
classical gravity spacetime of AdS, to be translated into information about 
dual defects probing properties of the conformal field theory. This has led to 
the determination of quantities like the energy loss of a moving defect, and 
the screening length between two oppositely charged defects. A well known 
subclass of the correspondence is that a plasma state of the CFT is dual to 
an asymptotically AdS space with a Schwarzschild black hole. 
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The current paper will examine the classical stability of some of the sta- 
tionary moving solutions found in [1] and will focus solely on the case of 
membranes in AdSr. As discussed in more detail there, this is an interest- 
ing case of the correspondence for various reasons. The dual six dimensional 
CFT is poorly understood, not having a Lagrangian description. The under- 
lying microscopic quantum gravity is M-Theory, allowing work in a less well 
tested case of the duality. The extended nature of the defect lends hope that 
it might explore different behavior of the CFT plasma than that probed by 
the zero dimensional endpoints of strings. The stability of these defects is 
important in this context; turbulent flow of the CFT plasma may manifest 
as instability of the probe. 

The paper is organized as follows: Section [2] recalls the infinitely extended 
moving membrane of [1] and its world-volume metric, and examines pertur- 
bations around this stationary solution, including the numerical calculation 
of the quasinormal modes (after a mathematical detour in Appendix [A] into 
multi-dimensional continued fractions); Section [s] discusses the special ana- 
lytic massless mode of the membrane corresponding to the Goldstone mode 
due to breaking of translation invariance; Section |4] recalls the compact mov- 
ing membrane of [1] and its world-volume metric, and examines the quasinor- 
mal modes of perturbations around this stationary solution; Lastly, Section 
[5] concludes with a discussion of the results. 

2 Infinitely extended moving membrane 

2.1 Background metric and equations of motion 

We will consider AdSy-Schwarzschild spacetinie with metric, G^^: 

d.s^ = — (-(l-u^)dt'+ ^ . (1) 

This metric is related to the Anti-dc Sitter black hole in the more familiar 
Poincare coordinates via u — . It has the horizon at u 1, and the 
boundary of the AdS space with radius of curvature i? is at u = 0. 

The equations of motion for the membrane are determined by the deter- 
minant of its induced metric, Qah = G^v^^a^\i where X>^[a) is the world- 
volume map into AdSr, and i(a\h denotes partial differentiation. From the 
Nambu-Goto-like action, S = —Tq J Y[i dai^J—g, where Tg is the membrane's 
tension, we obtain the equation of motion: 

(V^^?'^'G,,,x-,) ^ = . (2) 

2.2 Uniformly moving solution 

For a membrane extended along one of the transverse directions, y, and having 
a non-trivial profile in only one transverse direction, cc, we can parametrize 
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Figure 1: The x profile of tlie uniformly moving extended sheet as a function of 
the bulk radius u. 



the world-volume map as X^(a) = (t, u, y, x{t, u,y), z — zq) where z are the 
transverse directions upon which the map does not depend. For the ansatz 
of a stationary solution with no transverse dependence, x ~ vt + x{u), the 
equation of motion reduces to: 

x^ - ^^"^ f ] (3) 

where C is a constant that is determined by demanding that x\ and —g 
remain positive on the whole world- volume. This requires — 7^ — 1, 
where 7 is the Lorentz factor related to the "velocity" v. This reduces ([s]) to: 

V u 

21 — u-' 

which has the solution: 

Xy{u) = ±^ ^log I^a/i + U + v?^ - log(l - u) - v^arctan ^" ^^ , 

_(5) 

and is shown in Figure [T] The sign is fixed to be the lower one by examining 
the loss of energy into the black hole as the membrane moves along with the 
constant velocity v. 

2.3 World- volume metric and horizon 

In the original coordinates the world-volume metric gab is non-diagonal. By 
reparametrizing : 
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where 



5*11 



9tt (-l + r3) 
we have the non-zero components of the metric: 

i?2 



R2 



i?2 



9rr 



(6) 



In this form it is apparent that there is a horizon on the world- volume at the 
radial coordinate Uc = 7"^^^, where r = 1 and Qt-t and Qrr change sign. 



2.4 Expansion about stationary solution 

We will now consider perturbations to the above stationary solution. Again, 
the equations of motion will be determined by an action proportional to the 
square root of the determinant of the full induced metric, denoted hat- For a 
world- volume map X^(a) = X!^{a) + e''(a), with Xlj} the stationary solution 
above, and taking the perturbation and its derivatives to be small, we can 
expand the action of the perturbed solution to second order: 



+ 



d (d 



1 + 5^8? \ dXK \ dX% 



(7) 



The linear term can be rewritten as: 

'd_ 

) ^ 



dX. 



which vanishes after an integration by parts and use of the equation of motion 
for X^. 

The quadratic term can be rewritten as: 



d (d 



1 + 8^81 \ dX\ V 



Xi 



1 + si^S" V^"^ +GnsGi^-yXy,c^^^d\9 9 -9 9 -9 9 )) 



2.4.1 Transverse perturbations 

For a perturbation only in one direction transverse to the direction of motion 
we have e** = eS^, and the action for the perturbation becomes: 



1 + 



:^V-9Gzz9 



ab 
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Parametrizing the stationary world-volume with (r, r, w) gives a diago- 
nal g""^ and simplifies the covariant equation of motion for the transverse 
perturbation: 

(eb\/^G,,5'''),„ = . (8) 

For transverse perturbations of the form of a plane wave in the w direction, 
e(r, r, w) = g^T-^e'C^™^"'^', the equation of motion becomes the second order 
ODE: 



1 — r r ~ <j) r — 4 

'^"■^ _ _ ^3)^2^ = , (9) 



4r(l - r3)2 

where </> — e^'^^l'^ is a cubed root of unity. 

This equation is seen to have four singular points, tq, with the following 
respective characteristic exponents, A: 
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X 


0,2 


±iuj/6 


±10;/ (6(/)) 


±ia;/(6<?!>2) 



2.4.2 Boundary conditions 

In order to solve, or at least numerically find the quasinormal modes for 
the above perturbation we need to understand where and why to impose 
boundary conditions. As discussed above, from equation ^ for the world- 
volume metric, we can deduce the existence of a horizon on the world- volume 
of the stationary solution at the coordinate r — I. For densities and currents 
living on the world-volume (momenta etc.) r = 1 is a horizon and appropriate 
boundary conditions must be applied at this point, as discussed in [51 [3]. 

Near this world-volume horizon the local behavior of the perturbations 
is determined by the characteristic exponents at r = 1: A = ±iio/6. Their 
complex nature implies wave behavior in the coordinates r and t. At the 
horizon r = 1 we should have the boundary condition that we keep only 
those waves that are infalling, i.e. traveling to r > 1. A careful analysis is 
required because r diverges at ?' = 1, as we will now show. 

We are considering local solutions near the world-volume horizon of the 
form: 

e(T, r, w) = {l- r)±*"/6^(r)e'('="'-'^^) , (10) 

where S{r) is regular at the horizon. Near the world- volume horizon, for 
r < 1, this becomes: 

e « exp (— (=F log(l — r)/6 + t)) 

Recalling the definition of r we see that its near horizon behavior is: 



log(l — r) 1 — r 



r 



i2\ 
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This gives the near horizon behavior of the solution: 



We see that, for the top sign, the fogarithmic term is dominant, and following 
a constant phase, as t increases we require log(l — ?■) to increase, that is to 
move towards smaller r, away from the horizon at r = 1. This means the 
top sign corresponds to outgoing boundary conditions at the world-volume 
horizon. For the lower sign, the logarithmic term disappears and as t increases 
a constant phase requires (1 — r) to decrease, that is for r to get closer to 
the horizon at r = 1. Therefore the bottom sign corresponds to infalling 
boundary conditions at the world-volume horizon, and on physical grounds 
we will use the characteristic exponent —iuj/6 as a boundary condition for 
the perturbations. 

Making one last change of dependent function: 

e(r) = (1 - r)-*"/6(r - /'■^'I'^r - 

and defining the variable x = 1 — r, we obtain the equation of motion for 
r]{x): 

T]{x){k'^ - LO^) + ri'{x)A{tuj - 3 + 2(3 - tuj)x + {luj - 6)x^ + 2x^) 

+ T]" {x)4:{~3x + 6x^ ~ Ax^ + x^) ^ . (11) 

The above discussion of boundary conditions required rj to be regular at 
the horizon, so expanding in a local series solution about the horizon at a; = 0, 
r]{x) = X^^o''"'^"' ^® obtain the recursion relation for the coefficients c„, 
with n > 3: 

Cn-i {roSl + {n - i) {si + {n - i - l)U+i)) 

n(so + (n-l)i^^ ' ^ ' 



where rj, Sj, and ti denote the coefficients of x^ in equation ( 11 ) of ry, rj' , and 
rj", respectively. That is Tq = fc^ — w^, Si = 8(3 — iw), ^4 = 4, etc. The first 
three terms of the recursion relation are given by the initial conditions: 



U! 
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Co =1, Ci = — — , C2 



4(«w-3)' 32{ioj - 3)iiuj ~ 6) 

where cq is simply chosen as the normalization of the solution to our linear 
ODE. 

At the boundary, x — 1, the desired boundary condition is that rj vanishes. 
Demanding this and the ingoing wave condition at the horizon leads to a 
discrete set of possible w, the quasinormal frequencies. 
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2.4.3 Longitudinal perturbations 



For a perturbation in the direction of motion x, — eSf^, it can be shown 
that the action for the perturbation, equation Q, is simply 7^ times the 
action of the transverse perturbation. This was also the case for the string in 
AdSs, as seen in [3], and in fact is true for all the "analytically special" cases 
defined in [1] : those for which a defect of world- volume dimension n + 2 is in 
an AdS space of dimension d + 1 satisfy d = 2{n + 2) . 

Having an identical action up to scaling, for the longitudinal perturbations 
the equations of motion, and hence quasinormal frequencies, will be the same 
as for the transverse case. Properties dealing directly with the action, such 
as energy and momentum currents, will only need to be rescaled by powers 
of 7. The existence of the factor of 7^ can be understood in terms of Lorentz 
transformations. Going from the rest frame of the moving defect to the 
rest frame of the CFT plasma will lead to different transformations for the 
two point function of the longitudinal and transverse fluctuations. These 
correlators can be derived from the gravitational action under the duality, and 
it is seen that the 7^ factor is such to reproduce the correct transformation. 

That the factor of 7^ is the only difference between the two actions is 
a much less trivial statement. On the field theory side this is equivalent to 
saying that in the rest frame of the moving defect the two point functions 
for transverse and longitudinal fluctuations are equal. As pointed out in [5] 
for the case of a quark, this means that the fluctuations of the defect are 
spherically distributed in this frame, despite the effectively moving plasma 
background. This is only seen to hold for conformal theories, and is expected 
to be related to the fact that the drag coefficient is velocity independent in 
CFTs [Sj , although no general understanding from the field theory side exists 
at this time. 



2.5 Quasinormal frequencies 

Using the above equation of motion and boundary conditions, we can nu- 
merically determine the quasinormal frequencies in a variety of ways. The 
most straight forward is to numerically solve the equation of motion with the 
infalling boundary conditions at the world-volume horizon, for various values 
of k and w, and then determine for which values of Lu{k) are the solutions 
normalizable at the AdS boundary. This method is known to be limited in 
accuracy, especially at large uj and k. A similar approach, with similar prob- 
lems, is to truncate the power series solution for 77 at a finite number of terms. 



and using the knowledge of the coefficients from the recursion relations (12 1, 
we can also calculate for which values of uj{k) the solutions are normalizable 
at the AdS boundary. From past experience with the case of the trailing 
string in AdSs a more promising approach is to use continued fractions, and 
their relation to solutions of recursion relations. 
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2.5.1 Continued fractions 



An elegant and accurate method for obtaining the numerical values of quasi- 
normal frequencies implements the use of continued fractions, as discussed 
in [J. There, a criteria was used for which the local series solution had an 
increased radius of convergence, due to the existence of a minimal solution to 
( [T2| . This criteria was provided by Pincherle's Theorem which states that a 
minimal solution exists if and only if a certain continued fraction converges. 
By determining for which values of lj the fraction converges, the quasinormal 
frequencies are found. 



The recursion relations (12 1, unlike those of 0, are four term and there- 
fore Pincherle's Theorem, which only relates to three term relations, is in- 
applicable. Fortunately, due to [8] , a generalization to m -I- 2-term recursion 
relations exist, relating the existence of minimal solutions to the convergence 
of m-dimensional continued fractions. The pertinent results are discussed in 
the Appendix. 

In general, the continued fraction technique takes advantage of the fact 
that the characteristic exponents at the boundary differ by an integer, and 
hence the solution that behaves like contains logarithmic terms. Pincherle's 
theorem, and its generalization, gives us a criteria assuring that the recursion 



relation ( 12 ) has a minimal solution, which in turn assures that the solu- 
tion has an increased radius of convergence, and hence is analytic at r = 0. 
Therefore, for the values of uj{k) such that a minimal solution to the recur- 
sion relation exists, the corresponding solution to the equation of motion is 
analytic at the boundary, and hence must behave like r^; that is these uj 
correspond to Dirichlet solutions at the boundary with infalling boundary 
conditions at the horizon and hence are the quasinormal frequencies. 

The first few quasinormal frequencies for the infinitely extending mem- 
brane in AdS? with fc = are listed in the table below. From the differential 
equation it can analytically be shown that they come in pairs with real parts 
differing by a sign. The quoted significant figures are those which agree after 
running the continued fraction algorithm seventy and ninety steps (the a in 
/(°) in the Appendix); except the highest listed mode, for which the working 
precision of the numerics could not be made accurate enough after ninety 
steps. 



Re[Lo] 


Im[uj\ 


±3.8365833596 


-1.9990040306 


±6.51824614 


-3.52723336 


±9.1471822 


-5.036024 


±11.76121 


-6.5403 


±14.369 


-8.043 


±16.98 


-9.54 



Figure ([2| plots the quasinormal frequencies at fc = for various steps 
in the numeric iteration, showing the convergence behavior, which for the 
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Figure 2: The quasinormal frequencies in the complex oj plane, as determined 
numerically by the continued fraction method. Convergence is shown via plotting 
the frequencies as determined at differing number of steps in the iterative process. 
Shown are: Circles- 30 steps; Squares- 40 steps; Diamonds- 50 steps; Triangles- 
70 steps; Inverted Triangles- 90 steps. 



lower modes is much more accurate than the size of the dots. We note that 
all frequencies have a negative imaginary part, which indicate a decay of the 
perturbations and hence stability. The mode at the origin is simply the trivial 
uniform shift of the stationary solution x^. The structure in the complex 
plane of a near linear, equally spaced distribution is familiar from the case of 
perturbations in AdSs 

Figure ([3| shows the dependence of the lowest mode on the transverse 
wave number fc, for < k < Q. This behavior is also qualitatively similar to 
previous studies. 



3 The Goldstone mode 
3.1 An analytic solution 

Recalling the discussion above concerning boundary conditions at the world- 
volume horizon we are led to examine the equation of motion for 5(r), the 
regular part of the perturbation that has the infalling wave condition peeled 
off. Making the ansatz k = ±w for the dispersion relation this equation of 
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Figure 3: The transverse momentum dependence of the lowest quasinormal mode 
in the complex uj plane. As k increases from to 6 the mode moves up and to the 
right, although appears to stay bounded away from the real axis. 



motion has the analytic solution: 
S{r) =Ci exp 



— f — VSarctan ( 
6 V V 73 



C2 exp 



~6 



arctan 



1 + 2r 

~7r 



+ log 

log 



(1 - r)' 



(13) 



We note that this gives an e that is Neumann at the AdS boundary, r = 0. 
Examining the near horizon behavior we see that this solution behaves as: 

S{r « 1) = 3T5el^Ci + 3T^e^C2(l - r)"^ + O ((1 - rf) 

The second term interferes with the desired boundary conditions at the hori- 
zon and leads to the unphysical case of outgoing waves, we therefore require 
C2 = 0. Normalization then fixes C'l. Converting to our original variables, 
we can write the full perturbation as: 



6(a) = 
Ci cos 



-2/3 



±72/ ^ t ^ 7; ( arctan 



l + 2u 



log 



1 - u 



Vl + M + 



(14) 



This is plotted in Figure [l] for ?; = and ujj = 5. We see that in 
these field theory coordinates the perturbation has the effective frequency 
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t 



Figure 4: The analytic solution to the perturbation around the uniformly moving 
membrane. Plotted for y = and 077"^/^ = 5, with t E (0, 5) and u E (0, 1). The 
boundary is in the foreground at li = 0; from there the disturbance propagates 
causally down the membrane and bunches up on the bulk horizon at u = 1 where 
the radial wave number becomes infinite. 

oJeff = 7"^^"^^ and transverse wave number k^ff = ^^^^uj = "fuieff- From 
the boundary the perturbations propagate down the membrane, and bunch 
up at the bulk horizon u ^ 1, this being due to the increasing wave number 
effect of the logarithmic term log(l — u). 

3.2 Interpretation 

The dispersion relation in the field theory coordinates is easily interpreted. 
For a defect at rest in the plasma the dispersion relation is just that for a 
massless mode, i-^eff — ^eff- We can understand this as a Goldstone mode: 
the line defect in the field theory breaks the full translation invariance down 
to only along the y axis, the direction of its extent. This analytic solution 
then corresponds to the Goldstone modes for the four broken generators. For 
the moving defect a similar interpretation in terms of symmetry breaking 
exists, although the lack of Lorentz invariance due to the preferred rest frame 
of the plasma makes the dispersion relation less obvious. 



12 



4 Compact moving membrane 



4.1 Background solution 

We recall from [1] the solutions describing a membrane which is a circle on 
the boundary and moves with a constant velocity in a transverse direction, 
perpendicular to the enclosed disk. With this set up we can make the ansatz 
X = vt + x{u), and r = r{u) for the profile in the direction of motion, and the 
radius of the defect in the bulk, respectively. Written in traditional Poincare 
coordinates, u — z'^, the equation of motion for x{z) is very similar to that 
of section [2j and a constant of first integration allows us to write: 

^'"^ \l-zy (1 - z^)r^ - K^v^z<^ ■ ^ ' 

Here too the constant K is fixed by demanding x\ and —g are everywhere 
positive on the world-volume. This requires K — ±rc7, where = r^z^ = 
7^^/'^) is the radius of the defect at the critical bulk coordinate, Zc, where 
there is seen to be a world-volume horizon, as above. The sign is determined 
to be the bottom one, on the physical grounds that energy should be lost into 
the black hole, and not emerge. 

We can now substitute this into the equation of motion for r and obtain 
a second order nonlinear differential equation, which we numerically solved 
in [T] for different values of the velocity, u, and the critical radius, Tc- For 
Tc = 0, the profile must truncate and cap off smoothly at some bulk radial 
coordinate z < Zc, as seen by examining the energy- momentum flux down 
the world- volume. This gives the profile termed a "bowl" in 1 . For Tc > 
the defect stretches all the way to the bulk horizon, 2 = 1; solutions dubbed 
"tubes." For a range of parameters these tubes have unavoidable conical 
singularities where r — 0, see the discussion of [T]. 

4.2 World- volume metric and horizon 

For the reparametrization similar to above: 

p ^271/3, r^7-2/3(t + /(2)), 
we obtain a diagonal world-volume metric with; 

pj ^ ^"^ \pj (1 - 7- - ^c'"V 

The horizon at p = 1 is clear in these coordinates. As discussed above, the ex- 
istence of a world- volume horizon will be used to apply physically acceptable 
boundary conditions to perturbations about the above numerical solutions. 
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4.3 Perturbations 



4.3.1 Equation of motion 

We will consider perturbations about the above solutions for the radial em- 
bedding function; perturbations of the form e^{a) = e{a)6!^. Considering e 
and its derivatives small, we can use the expansion of the action to second or- 
der, equation ([t]), to obtain the equation of motion outside the world- volume 
horizon: 

e"{p) + P^e'ip) + (~m''- co^^-) e{p) = , (16) 
P \ P Pj 

where we have used the parametrization and ansatz e(a) — exp[zm0— ia;r]e(/9), 
and: 

r2 + rlv'^p'^r'^^ I (1 - p6){{l - -f-^p(i)r^ ~ rlv"^ p^) 




(1 — 7 '^p^)r'^ — r'^v^p^ ' 
r^(l + (l-7'V)rfJ 

(1 - /96)((1 _ ^-2^6)^2 _ ^2^2^6) 

4.3.2 Boundary conditions and quasinormal modes 



(17) 



As discussed above for the case of the infinitely extended membrane, the 
existence of a world-volume horizon gives us a location and rationale for the 
imposition of boundary conditions. At this horizon we should have infalling 
boundary conditions on our perturbations. We can understand how this will 
manifest by examining the equation of motion ( |16[ ) near this singular point. 
A power series for r near p = 1 was determined in [T] in order to implement 
the numerics finding the background solutions. This analytic structure can 
be used to determine the characteristic exponent near the horizon. For e « 
(1 - p)^ it is found that: 



The complex nature of the exponent leads to wave behavior in the p — t 
coordinates, and demanding infalling boundary conditions picks the lower 
sign. 

We now write e{p) = (1 — p)^^^^'q[p) where 77 is regular at the horizon. 
The known power series approximation for r now allows the determination 
of the power series behavior of 77 near the horizon. As in [1], this gives 
us enough boundary conditions at the horizon to numerically integrate the 
equation of motion. To determine the quasinormal frequencies we examine 
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-0.5 - . 

-1.0 - ' . 

-1.5 - 
• • • • -2.of • 

Figure 5: The lowest quasinormal frequency of the tube solutions, with 7 = 6 
and m = 0. From top to bottom Tc increases from .002 to 1.1, approximately in 
logarithmically spaced steps. 



the boundary behavior of rj and numerically determine for which frequencies 
uj(m) the solutions have the desired vanishing Dirichlet values. Unlike in 
section [2] the numerical nature of our background solution precludes the use 
of any continued fraction method. The lowest lying quasinormal frequencies 
are shown in Figure [5] for 7 = 6, m = 0, and various values of rc', Figure 
[6] shows the m = 1 case. The numerics much beyond these lowest lying 
modes is severely limited in accuracy. As expected, in the large Vc limit, 
these frequencies approach that of the infinitely extended sheet of Section |2] 
Of consequence is the non-positive nature of the imaginary part of the 
quasinormal frequencies. This indicates that the perturbations decay in time, 
and that the background tube solutions are classically stable. This comes 
as somewhat of a surprise, as the situation is reminiscent of the Gregory- 
Laflamme [10] situation, where a classical instability is seen for extended 
objects (in their case black strings) which have a proper radius smaller than 
the radius of curvature of the AdS space. In our coordinates this corresponds 
to the radius of the tube of order the bulk radial coordinate z; in terms of the 
radius at the world-volume horizon this translates to rc = Zc, or — 6^^/"^ ~ 
0.55 for the case shown in the figures. Nothing special appears to happen in 
this regime; although the quasinormal frequencies approach the real axis as 
Tc decreases, indicating an increase in the characteristic damping time, they 
do not cross the axis and lead to unstable growing modes. 
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Figure 6: The lowest quasinormal frequency of the tube solutions, with 7 = 6 and 
m = 1. Again Tc increases from top to bottom. 



5 Conclusion 

The major purpose of this paper was the determination of the stability of 
moving membrane configurations in Ad.S7. For the infinitely extended "sheet" 
of Section [2] we used a continued fraction method to numerically find that 
the quasinormal frequencies of the membrane all have a negative imaginary 
part. Similar conclusions were drawn for the compact "tube" solutions of 
Section |4] using numerical integration. Perturbations to the compact "bowl" 
solutions are trivially stable in this sense: having Dirichlet conditions at the 
boundary and regularity conditions at the bottom of the bowl, instead of 
the inf ailing wave conditions for the tube, they will have fully real normal 
modes, as was confirmed with numerics. Thus we deduce the stability of 
all the classical solutions found in [T]. This was somewhat unexpected for 
the case of the tubes, where a Gregory-Laflamme instability was expected. 
Additionally, the hope of observing hints of turbulent behavior of the plasma 
fiow through an instability of the probes is not borne out. 

Two unexpected finds surfaced during this work. The first is the determi- 
nation of the massless Goldstone modes of Section [3] These can be expected 
to exist on symmetry breaking grounds, but to have found an analytic solution 
for their behavior is a lucky event. A fuller understanding of their properties 
would be desirable. The other interesting find was the relation between the 



actions for the longitudinal and transverse perturbations, discussed in 2.4.3 
Thus for the analytically special cases of d = 2(n + 2), we have shown that 
in its rest frame an rt-spatial dimensional defect in the d-dimensional CFT 
has a symmetric distribution of its perturbations, despite the presence of an 
effectively moving plasma. 
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The last interesting development was the use of higher dimensional con- 
tinued fractions to calculate quasinormal modes. This approach should be 
applicable to any set of /c-term recursion relations, and not just the three-term 
case that has been most exploited in the literature via Pincherle's Theorem. 
We believe that this is the first time Parusnikov's generalization has been 
used in this context, and expect it to have much a greater application. 
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A A generalization of Pincherle's theorem 

In [5] the author generalizes Pincherle's theorem to relate TO-|-2-term recursion 
relations and m-dimensional continued fractions; that is vectorial continued 
fractions. This necessitates the use of projective vectors to accomplish "di- 
vision by a vector." The needed results, along with required definitions, are 
discussed in the following. 

Parusnikov's Theorem 2 assures us the following correspondence: 



Theorem 2- 



m+l 



iii) (to -|- 2)-term recursion relation: qn = Pm+2-j.nQn-j 



i) interruption-free TO-dimensional continued fraction: 

ao 



/o - foo + 



bi 



(18) 



where p„ = (a„ ^, a„ ^&i^„, • • • ,a„^6m,ri); "interruption-free" means the n-th 
remainder, a„/(6„_|_i -I- • • • ), is non-zero; and throughout the notation Vi,j will 
denote the i-th component of vector Vj. 

The construction of higher dimensional continued fractions involves pro- 
jective vectors, and linear algebra. Given the vector p„ defined above, a 
simple iterative algorithm for numerically calculating continued fraction /o 
follows: Define: 



P2,n 



\ 



V 



1 P: 



17 



/ 1 
1 



\ 





... 1 
\ 1 ... / 

These matrices act on the m + l-dimensional projective space, denoted 
(/i : ■ • • : /m+i) with fi € C, and we will be interested in the subspace of 
m + l-dimensional vectors whose last component does not vanish: 



Lastly we need the map: 



fm+l ^ 0} 



>?Tl+l 



We then have: 



frn+l) 



hn+1 



hn+1 



fo = lim a ^QiCm+i 



where e 



m+l 



(0 



0:1) 



This theorem provides a mapping between recursion relations, such as 
those arising from differential equations (12), and multidimensional contin- 
ued fractions. The next theorem provides a powerful result concerning the 
convergence of a continued fraction and the existence of minimal solutions of 
the recursion relations. 

Theorem 6— 



An m-dimensional continued fraction fo converges. 



There exists a minimal subspace of the space of sequences satisfying 



the corresponding recursion relations given in (18 1. 
Where we recall that a sequence, /i„, is minimal if: 

Hm ^ = 

n^oo g„ 

for all other sequences (/„. 

Theorems 7 and 8 then give a recursion relation satisfied by the minimal 
solution, allowing an explicit construction of a basis. 

Theorem 7/8- 

The minimal basis obeys: 

, _ / ~fj,n+l \ , f -fm+l,n+l , , 



m,n+l / V frn.,n-\-l 
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The notation here is a bit opaque, and the interested reader is referred to 
the original work [8]; the apphcable statement to the case at hand fohows. 
Additionally, this formula corrects a typo in the original, in which the roles 
of fm+i,n+i and /m,n+i are reversed. This is required to be consistent with 
the line preceding Equation (23) in [5], as well as with the check of part 
(B) of Pincherle's theorem following Theorem 8. As a final assurance, this 
correction is needed for the continued fraction results to agree with other 
numerical methods. 

For the case at hand we can summarize these results into a compact 
formula used to numerically calculate the quasinormal frequencies. Given 



the recursion relations (12) and the initial conditions c; for I = 0, 1,2, then 
for: 

Parusnikov's Theorem 7 gives: 

C2= lim f(^)ci+|-^)co) , (19) 




for the minimal solutions, subject to the above initial conditions. 

We are interested in the minimal solutions to our equation of motion 
because these solutions are such that their radius of convergence is increased, 
and hence they are analytic at the singular point of the boundary, r = 0. 
Recalling that the characteristic exponents at the boundary are (0,2), we 
know near boundary solutions will behave as: 

yi{r) ^ + 0{r^) 

y2 ^ I + bir + h\og{r)yi{r) + 0{r^) 

where bi{u),k) and /i(aj, fe) are determined by the recursion relations (12), and 
the crucial logarithmic term arises due to the characteristic exponents differ- 
ing by an integer. Thus we see that the solution with non-Dirichlet bound- 
ary conditions at r = is also the solution with a non-analytic logarithmic 
term. Therefore, the minimal solution, being analytic at the boundary, is 
also Dirichlet, and given its infalling behavior at the horizon corresponds to 
a quasinormal mode. This generalization of Pincherle's theorem gives us a 



criteria for determining those uj such that a minimal solution to ( 12 ) exists. 



i.e. quasinormal frequencies are those for wh ich ( 19 1 holds. Numerically, we 



calculate /^"^ for large a and then use (19) to calculate w. Repeating for 
larger a gives us an understanding of the convergence properties and allows 
the approximation of the quasinormal frequencies to a desired accuracy. 
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